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We present evidence for collective action of Cr^^ ion impurities in a highly concentrated ruby crystal 
coupled to microwave Whispering Callery Modes (WCMs). The cylindrical geometry of the crystal 
allows for the creation of superradiant, or “spin-mode” doublets, with spatial structure similar to 
that of WCMs. The formation of these spin patterns allows us to observe directly different selection 
rules namely wavenumber and azimuthal phase matching. The demonstration is made via an avoided 
level crossing between spin and photon mode doublets as well as absence of coupling between spin 
modes of different wavenumbers. The effect is observable due to strong spin-photon coupling (67 
MHz) exceeding both spin ensemble and cavity losses as well as the photon doublet splitting. We 
demonstrate that a four harmonic oscillator model not only with coupling between photon resonances 
( 0.43 MHz) but also with spin doublet (73 MHz) is necessary to accurately describe these results. 


S uperradiance is an important phenomenon in quan¬ 
tum optics as often the sample under study features 
separation distances small compared to the wavelength of 
exciting radiation, A. There is renewed interest in these 
systems for quantum information sciences and to attain 
new insights into QED and its applications [1, 2]. Eor 
example, superradiant effects will need to be considered 
when constructing an optical-microwave interface [3-5] 
or a quantum memory [6, 7] using spin ensembles. 

Superradiance was initially defined in 1954 by Robert 
Dicke as the cooperative, spontaneous emission of pho¬ 
tons from a collection of atoms [8]. Superradiance 
was first observed experimentally in 1973 in the opti¬ 
cal regime in HE Gas [10]. It has since been observed in 
other ultacold atomic gases [11-14], organic semiconduc¬ 
tors [15, 16], polymer thin films [17], numerous crystalline 
systems [18-20], and in artificial atoms [21-23]. Here, we 
report the observation of superradiance in the microwave 
regime in a highly doped ruby sample, with relatively 
high concentrations of Gr^+ ions replacing Al^+ ions in 
the crystal lattice. 

In free space, when N atoms are close together com¬ 
pared with A, they act like one big atom and decay collec¬ 
tively, in phase with one another. As a result, the atoms 
radiate their energy N times faster than for incoherent 
emission. A direct result is the inherent directionality as¬ 
sociated with the emitted radiation; the emitted photons 
travel in the same direction as the exciting photons. This 
directionality is a result of the timing of the excitations; 
the atoms at the “front” of the sample are excited first, 
and those at the back, last, leading to the excitations ap¬ 
pearing as spatial phase factors [9]. Superradiance is a 
consequence of extra coherence in the system, which can 
be observed in additional ways on top of an increased 
emission rate. 

To observe coherent effects originating from collective 
action in the microwave regime is sufficiently more 


challenging than the optical regime. This is due to 
the relatively weak strength of field-matter interactions 
via magnetic fields as compared to electric fields [24]. 
When the emitters couple to a resonant cavity mode, 
their separation becomes irrelevant. The coherence 
between separate spins is generated by their interaction 
with a common mode, which occupies space over the 
entire cavity volume. The strength of this interaction 
is determined by the light-matter coupling constant, 
g, which is proportional to the concentration of the 
emitters [25]. 


Unlike Eabry-Perot cavities, ideal whispering gallery 
mode (WGM) resonators have rotational symmetry. This 
fact dictates that if a mode field distribution has solutions 
of the system eigenvalue problem, any of its rotations 
around the cylinder axis will also be a solution. Each of 
these solutions could be represented as a linear combina¬ 
tion of only two orthogonal solutions. In actual WGM 
cavities, this symmetry is lifted by a number of imper¬ 
fections that we further collectively call back-scatterers. 
These back-scatterers introduce a coupling between the 
two particular orthogonal solutions that depend on the 
back-scatterer details. A WGM will therefore manifest 
as two orthogonal modes (or doublets) with a difference 
of Sine and Gosine in the mode’s azimuthal dependence 
in its analytical expression (i.e. a difference of 7r/2 in az¬ 
imuthal phase) [26], henceforth referred to as the “ 5 ” and 
“c” modes. This manifests as a splitting of a single reso¬ 
nant peak into two resonant peaks by a distance equal to 
two times the coupling value, n. In sapphire crystals, the 
losses of such WGMs are so low, that the bandwidth of 
these modes is generally less than 2/^ hence the doublet 
resonance can be resolved. 

The Hamiltonian describing such a WGM doublet res¬ 
onance is 
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Here ujk is the angular frequency of a WGM with 
wavenumber and a|.^, a/^^s, a|. ^ and ak^c are the 
bosonic raising and lowering operators of the “s” and 
“c” doublet constituents of this WGM, respectively. The 
first term in equation (1) represents both modes as sim¬ 
ple harmonic oscillators (SHOs), while the second term 
represents the coupling between them, which produces 
the mode-splitting and doublet appearance. 

A crystal containing dilute concentrations of param¬ 
agnetic ion impurities will demonstrate an absorption of 
energy from these WGMs into the spin angular momen¬ 
tum of the ion’s valence electrons if the frequency of the 
latter transition is tuned (via the Zeeman effect) to be 
coincident with that of the former. Only WGMs with 
magnetic field components perpendicular to the applied 
DG magnetic field will interact in this fashion. This 
limits the discussion to WGMs that are polarised with 
a field distribution (“WGH” modes), since 

the applied magnetic field in the described case is aligned 
with the z-axis of the crystal. 

In general, the collective electron spin resonance (ESR) 
can be considered as an ensemble of independent, non¬ 
interacting two level systems (TLSs); and the crystal it¬ 
self as a paramagnetic material. In such a case, the ESR, 
WGM doublet and the interaction between the two can 
be described by the modified Tavis-Gummings Hamilto¬ 
nian: 


Htc = Ho+ 'E + 
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E E ) + 
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Here, uoi is the angular resonant frequency of the 
TLS transition at a particular B-field, and aj' and 
are its raising and lowering operators, gk^s {gk,c) is the 
coupling between the spin ensemble and the “ 5 ” (“c”) 
mode of the WGM with wavenumber k. One of these 
coupling terms will be set to zero because the ESR can 
only couple to one resonance of the doublet due to the 
spin conservation law [27]. Note that the choice of which 
mode will couple (“s” or “c”) is determined by the sign 
of the change in spin angular momentum of the TLS 
transition in question; Am = ±1. 



EIG. 1: (Golor online) Eour harmonic oscillator model 
depiction. 

A spin ensemble can be considered as a classical sys¬ 
tem of SHOs distributed over a large region of space. 
Eor densely packed ensembles interacting with a com¬ 
mon cavity mode, these SHOs can be phased relative to 
each other so that coherent radiation is obtained in a 
particular direction. This is referred to as superradiance, 
and occurs when a group of N emitters interact with a 
common light field in a collective and coherent fashion 
[ 8 ], 

The coherent radiation generated by excited atoms 
emitting photons, hereon in referred to as the “spin¬ 
mode” , is completely analogous to a photonic WGM, and 
therefore can exist as a doublet due to backscatterers, in 
exactly the same way. It will also display the same type 
of wavenumber orthogonality, doublet orthogonality and 
a coupling between the two doublet constituents. Spin 
doublet modes have been previously observed in ferro¬ 
magnetic YIG samples [28], but never before in a doped 
sapphire system. In such a scenario, the Hamiltonian 
describing the interaction between the photonic cavity 
WGM doublet and the spin-mode doublet would appear 
as 

H = Ho -^'^Ldk (^<^k,s^k,s + ^k,c^k,c^ 

k 
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where y represents the coupling between the two spin¬ 
mode doublet constituents; “ 5 ” and “c”. This Hamilto¬ 
nian is derived following the treatment of Dicke [8] when 
describing radiation from a gas of large extent. Erom eq. 
(2); a summation over all modes and TLSs, a transition is 
made to just the former. The selection rules of such a sys¬ 
tem [8] dictate that only modes with equal wavenumbers. 
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/c, may interact. In addition to this, the equivalent dou¬ 
blet orthogonality of the spin-modes and WGMs allow 
for both and spin-WGM interactions, 

but not This removes the requirement that one 

of the spin-mode couplings be set to zero. It is reason¬ 
able to assume that the coupling strengths of the two “s” 
modes will be equal to the two “c” modes, hence the use 
of a non-polarisation-specific coupling term, gk. As such, 
the allowed spin-WGM interactions are described by the 
third term in eq. (3). 

The second expression in (3) represents the spin-mode 
doublets as two SHOs, while the final term represents the 
coupling between them, resulting from imperfections in 
the crystal; a direct analogue of the last term in eq. (1). 

Equation (3) is represented diagrammatically in FIG. 
1 for a single value of k. It describes a scenario of 
four SHOs with the allowed linear couplings. This is 
distinctly different from the case described by eq. (2), 
which would exists as three SHOs [27]. 



FIG. 2: (Color online) Experimental setup. The crystal 
and superconducting magnet are cooled to 4 K in a 
cryogenic refrigerator, whilst the microwave coupling 
probes are aligned to excite WGH modes {Ez^ H(jy 
components). Two probes are used to view the relevant 
modes in transmission (S21) on a Vector Network 
Analyser (VNA) from which Q factors can be 
determined as well as the frequency shift caused by 
coupling to spins. 


The experimental set up is identical to that described by 
Farr et al. [25], however we examine WGMs closer to the 
zero-field splitting levels of the Cr^+ ensemble. The ori¬ 


entation of the crystal, microwave coupling probes and 
applied DC magnetic field is depicted in FIG. 2. Typi¬ 
cal ESR parameters for Cr^+ ions can be found in [29]. 
In this paper, we deal with the Am = ±1 transitions; 
|—3/2) ^ \ — l/2) and |3/2) ^ |l/2). Both these tran¬ 
sitions have a zero-field frequency of 11.447 GHz, and 
tune in opposite directions as B field is swept (Am = +1 
increases in frequency with an increase in B field, and 
vice versa) with df/dB = where is the Lande 

g factor and (3 is the Bohr magneton. 



Magnetic Field (T) 


FIG. 3: (Color online) Spectroscopy results showing the 
interaction of the generated and pump WGMs with the 
1+3/2) ^ 1+1/2) Cr^+ electron spin transition as 
magnetic field is swept. 

For example, FIG. 3 shows the Am = — 1 transition 
(in red), as it moves through five distinct WGMs (in 
black). Each of the data points that make up the black 
curves represent the position of a resonant peak within 
a single 8 MHz sweep centred around that particular 
frequency for that particular magnetic field value. The 
power incident on the crystal is Pine = —60 dBm, which 
corresponds to a photon occupation number on the 
order of 10^. Far from the intersection of the WGMs 
and the ESR transition, the black curves in FIG. 3 
represent the frequency location of the WGMs, however 
when the ESR is tuned such that a particular WGM 
is within its bandwidth, the black curves represent 
hybrid spin-WGMs, and an avoided level crossing (ALG) 
can be observed, as depicted in the inset figures of FIG. 3. 

Equation (2) predicts a gyrotropic response for the ALG 
of a WGM doublet and ESR, which may be modelled 
with great accuracy by three SHOs [27]. The ESR spec¬ 
troscopy results for the ruby crystal in question (FIG. 
4 and 5) clearly show an absence of this gyrotropic re¬ 
sponse. We observe that both components of the WGM 
doublet interact with the spin transition. 

As FIG. 4 demonstrates, there are four asymptotes to 
which the hybrid modes converge. The two horizontal 
asymptotes of FIG. 4 are a standard result of the WGM 
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FIG. 4: (Color online) Interaction of the 9.545 GHz 
WGM (WGH 7 ,i,i) and the |+3/2) ^ |+l/2) 
transition at H = 0.68 T. Asymptotes (dashed) and four 
SHO model fit (solid). 


in question (/ = 9.55 GHz, transverse magnetic mode 
with 7 azimuthal nodes in 180° and 1 radial and axial 
node - WGHyqq) existing as a doublet. The vertical 
asymptotes, which in fact depend on B-field just as the 
red curve in FIG. 3 (observable if the y-axis scale were 
broader), confirm the presence of a spin-mode doublet. 
Their separation is 2%. In the general paramagnetic case 
of equation (2), there would be only one vertical asymp¬ 
tote [27]. 

The presence of these four asymptotes requires a four 
SHO model (FIG. 1) to fit the experimental data. Using 
values of of ^ = 67 MHz, k = 0.43 MHz and y = 76 MHz, 
a fit is produced which is displayed in FIG. 4. To produce 
good agreement with the model, it is essential that the 
cross coupling terms gx be neglected, or at least be much 
smaller than the spin-mode couplings g - consistent with 
the allowed terms in the third expression in eq. (3). 

The requirement for spatial orthogonality of the 
spin-modes is again confirmed by FIG. 5. Here, we 
see the same type of ALC as in FIG. 4. However, we 
also observe the tail end of another ALC originating at 
a slightly higher frequency enter the frame. It is the 
hybrid mode of a higher frequency spin-mode and WGM 
of a different order. As predicted by the selection rules 
of such a system, these two doublets simply merge; there 
is no interaction, due to their different wavenumbers. 

As WGMs hybridise with a spin-mode (or paramagnetic 
spin ensembles, for that matter), not only is a frequency 
shift observable due to the altered magnetic susceptibil¬ 
ity of the resonant dielectric, but a change in the now 
hybrid mode Q factor becomes apparent. This is due to 
an additional loss mechanism introduced via the coupling 
to the spin-mode. 

The Q factor of the 9.5 GHz mode as B-field is tuned 
through the ESR centre is plotted in FIG. 6 for both 
constituents of the mode doublet. As the ESR becomes 
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EIG. 5: (Color online) Interaction of the 12.74 GHz 
WGM and the Cr^+ l~3/2) ^ | —1/2) transition. 



Magnetic Field (T) 

EIG. 6: (Color online)Q factor of the hybrid “ 5 ” and 
“c” modes at 9.54 GHz as B-field is swept. 

more closely tuned to the WGM frequency, the mode 
hybridises to a greater extent and Q factor drops. This 
is evidence that the losses experienced by spin-modes 
are greater than those of the photons. Q factors were 
measured by fitting a Eano resonance line shape [30] to 
the S21 data obtained from the VNA at discrete B-field 
values. 

The four HO model fit to the experimental data allows 
one to make an estimate of the concentration of Cr^+ ions 
in the ruby [25]. The concentration of spins participating 
in the WGM interaction, Uparti can be calculated using 
the value of the spin-mode coupling, g and the magnetic 
filling factor of the WGM in directions perpendicular to 
the applied DC magnetic field, Eor the 9.5 GHz WGM, 
^ = 0.877 and was calculated from finite element mod¬ 
elling. Tipart is calculated to be 8.34 x 10^^ ions/m^. 

At finite temperatures, in the continuously driven 
regime, the number of ions prepared in either the ground 
{N-) or excited states (A^+) of the relevant transition is a 
function of the thermal distribution of ions and the total 
number of impurity ions in the crystal, Nt- By equating 
^partV = A"-, where V is the total volume of the ruby 
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crystal, one can solve for Nt {Nt = 5.42 x 10^^ ions). 
Given that there are two Al^+ ions per unit cell of sap¬ 
phire that Cr^+ can potentially replace, and the volume 
of the unit cell is that of a trigonal crystal system, the 
total concentration of ion impurities can be calculated 
as Nt divided by the total number of potential lattice 
cites for Cr^+ ions to take. A concentration of approxi¬ 
mately 40 ppm Cr^+ is calculated. This agrees very well 
with previously measured values for this same crystal, re¬ 
ported as 34 ppm [25, 29], hence confirming the validity 
of the four SHO approximation used here, and ergo the 
conclusions that can be drawn from it. 

This concentration is approximately two orders of mag¬ 
nitude larger than the concentration of Fe^+ impurities 
in [27] (150 ppb), and the value of g is also approximately 
an order of magnitude larger. In addition, the losses as¬ 
sociated with the Cr^+ ESR {Aujspins = 9 MHz) are 
three times less than those associated with the Fe^+ case 
{Aujspins = 27 MHz) [27, 31]. In the present case, 
g > AuJspins, AuJwGM (AcJvVGM/27r = Q/fres = 1-6 
kHz), satisfying the conditions for strong coupling. It 
is due to this strong atom-field coupling that superradi¬ 
ance can occur as it is this interaction from which the 
collective action of the ensemble is derived, and large ion 
concentrations contribute to this (as g = gQ\fN). This 
explains why a superradiant ESR, and hence a four SHO 
model, is observed in the Cr^+ case and not in the pre¬ 
viously reported Ee^+ case [27, 31], which did not satisfy 
the conditions of strong coupling. 

The tell tale sign of strong coupling (splitting of the 
resonant cavity mode when the ESR is tuned) is however 
unobservable due to the loss of coupling between the 
microwave pump source and WGM when the two are 
tuned. Because cavity losses are orders of magnitude 
lower than those associated with the spin ensemble, 
and the coupling between the transmission line and 
the WGM, yd, is proportional to Q, when the ESR is 
tuned to the WGM frequency the extra dissipation 
which is introduced drastically reduces the ability of 
the transmission line to excite the WGM. Despite the 
fact that the sum of the spin and cavity mode losses 
is less than the spin-mode coupling, this change in 
external coupling to the input microwave probe results 
in the mode vanishing. This is why no resonant mode is 
observed near the ESR centre as shown in the insets in 
EIG. 3, in EIG. 4 about AB = 0, and hence why no Q 
factors can be derived in this region in EIG. 6. 

We have described a study regarding the collective inter¬ 
action of Ghromium impurity ions in ruby with crystal 
photonic WGMs, resulting in the creation of a spin¬ 
mode doublet. Theoretical predictions and experimental 
measurements reveal a set of selection rules that govern 
the interaction between the spin-modes and WGMs: 
wavenmuber and azimuthal phase matching. We observe 
an avoided level crossing between WGM and spin-mode 


doublets in a fashion that can only be described by a four 
harmonic oscillator model. These results may have im¬ 
portant implications for QED experiments dealing with 
strongly coupled light-matter interactions, as this new 
model must be used to describe the resulting phenomena. 
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